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PREFACE 


The purpose of this report is to present a technique of engineer- 
ing decision making. The technique will employ the engineer's expert 
opinions in the form of probabilities concerning unknown conditions, 
and his knowledge of costs or values if these possible conditions 
were to exist. 


Part I of this report deals with Posterior Analysis, which is 
selection of an action decision where all opinions and investigation 
results are available, and the probabilities of unknown conditions 
have been evaluated with this information. 


Part II deals with Pre=Posterior Analysis in which the engineer 
feels that he does not yet have enough information to select an action, 
and therefore he wishes to choose a level of investigation to provide 
results which will better enable him to evaluate the probabilities 
of unknown conditions. The first example in Part II (extension of 
the umbrella problem) serves to introduce the concepts of Pre-Posterior 
Analysis. This example is followed by a sophisticated problem, the 
solution of which illustrates a formal and somewhat general procedure, 
and thus is not essential in a first reading. 
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"Manuscript" Statement 


This edition is considered a rough draft and subject to correction 
before final incorporation into a Forest Service Technical Report. 


This edition will be revised about January 1975 for final incorpora- 
tion into the Forest Service technical report system. Interested 
parties should advise the Division of Engineering at the Forest 
Service National Headquarters of any changes they would like to 

see made in this document. Such changes will be considered during 
final revision, 


We believe that techniques discussed in this edition are accurate, 
reliable, and the best available for the use listed. Therefore, a 
field unit should not hesitate to use the techniques. 


We expect that changes will generally be in the area of clarification 
or in correcting obvious errors, rather than in changing the techniques. 
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DECISIONS CONCERNING THE SELECTION OF AN ACTION IN THE FACE OF 
UNCERTAINTY (POSTERIOR ANALYSTS) 


The necessity of making decisions in the face of uncertainty occurs 
with disturbing regularity in everyone's personal and professional life. 
The engineer, in particular, must make decisions which involve significant 
economic consequences with respect to project planning and selection of 
investigations related to future projects. He must employ all of his 
professional knowledge concerning techniques of forecasting, reconnais- 
sance, and exploration in order to reduce the uncertainties related to 
future requirements, and his knowledge must also provide workable esti- 
mates of the costs and benefits related to any proposed project. He 
must then employ all of this knowledge to select the best type or size 
of project. A pertinent example involves the choice of the reliability 
level of flood protection for a given area: the engineer must decide 
whether a 90 percent reliable project with cost of one million dollars 
is more advantageous than a 99 percent reliable project with cost of 
on million dollars, when the loss due to flooding may be one million 

ollars. 


The purpose of this report is to present a logical method of making 
consistent decisions. The method will employ all of the engineer's pro- 
fessional knowledge concerning the probability of future or unknown con- 
ditions, and information related to costs and benefits. This method is 
based on the concept of minimizing the expected loss, or maximizing the 
expected gain, and can be introduced by considering the following situa- 
tion. 


Let us suppose that we are invited to indulge in the following game 
of chance: a die is rolled and if ace occurs, we win one hundred dollars, 
and if not ace, the we lose one dollar. Very few, if any of us would 
refuse this fairly good chance of winning one hundred dollars because 
of the relatively small amount of possible loss. Our decision to accept 
this attractive gamble is based on a thinking process which is mathemat- 
ically equivalent to a computation of expected gain: defined as, 


(winnings) x (chance of success) = (loss) x (chance of failure). 
For example, in our game the expected gain is 
($100)(1/6) = ($1)(5/6) = $15.84 


It can be shown that a decision based on the maximum expected gain, 
or minimum expected loss, is economically sound. In engineering, economic 
decisions involve gambling in the face of uncertainty quite similar to 
our proposed game of chance, and therefore it is necessary to develop a 
minimum expected loss calculation procedure, involving engineering esti- 
mates of probabilities or chances and costs, for any given decision 
problem. This general decision procedure or method will be developed 
in terms of the following, very human, example in which the losses and 
probabilities are readily available: 
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A cosmopolitan gentleman has deemed it absolutely necessary 
to attend a theatrical performance. He must walk home from 
the theatre, and therefore consults the newspaper weather 
report for the possibility of rain at the end of the perfor- 
mance. The current method of weather forecasting furnishes 
the gentleman with a 40 percent chance for rain at the time 
that he must walk home. His decision in the face of this 
uncertainty of rain is whether or not to carry an umbrella. 
If he goes without the umbrella, and rain occurs, then he 
must pay two dollars to press his rainsoaked suit. If he 
elects to carry the umbrella, then he must pay a fifty-cent 
fee to check in the umbrella at the theatre cloakroom. 


I(a) Symbols and the Decision Tree 


In order to best present all of the information in this problem, 
and to provide a general calculation procedure in terms of this infor- 
mation we will define the following symbols and introduce the concept 
of the decision tree. 


List of Decision choices or Actions {1,2} 


1 


Take umbrella 
2 = Do not take umbrella 
List of Future Actual Conditions or Situations {A,B} 
A = rain at end of theatrical performance 
B = no rain at end of theatrical performance 
Probabilities of Future Actual Situations P(A), P(B) 
P(A) = probability of rain = 0.4 
P(B) = probability of no rain = 0.6 


Losses (L) corresponding to a selected action (1 or 2) and actual 
situation (A or B). 


L(1,A) = Loss if umbrella is taken (action 1) and rain occurs (situation A) 
= $0.50 cloakroom fee 


L(1,B) = Loss if umbrella is taken (action 1) and no rain occurs (situation B) 
= $0.50 cloakroom fee 
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L(2,A) = Loss if no umbrella is taken (action 2) and rain occurs (situation A) 
= $2.00 suit pressing charge 


L(2,B) = Loss if no umbrella is taken (action 2) and no rain occurs (situation B) 
= 0 
The above information can be presented in a diagram termed a decision 


tree, which allows an orderly presentation of the action choices 
and the loss consequences due to future uncertain actual conditions. 
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Here the gentleman, as the decision maker, stands at the trunk of 
the tree: the main branch system displays the actions which may be 
selected; and the secondary branch system displays the possible future © 
conditions along with their corresponding probabilities and loss conse- 
quences. This decision tree presentation allows the decision maker to 
make an orderly visual analysis of any of the different combinations 
of actions and actual future conditions. For example, if we follow 
along the action 1 branch, we have a complete display of both the 
probability of outcome A, P(A) and the loss consequences L(1,A) if 
action 1 is taken and A occurs. 


In actual engineering decision situations, where more than two 
actions and conditions may occur, the decision tree offers the best 
means of presenting all information to a group of decision makers. 
Each action, actual condition, probability, and loss consequence 
can be discussed individually and compared with the other information 
on separate tree branches, 


(b) Choice of Action based on Minimum Expected Loss 


In our simple example, the gentlemen will select the best action 
based on the criterion of minimum expected loss. Each action (J) has 
an expected, or average loss given by 


L (J) = Sum of the Products of the losses (consequences) which may 


occur if action J is chosen multiplied by the probabilities 
of their occurrence. 


L(J) = L(J,A)P(A) + L(J,B)P(B) 


For our particular problem 


qa 
uou 


1, Take umbrella, and 
2, Do not take umbrella, 


the expected losses are: 


L(1) = L(1,A)P(A) + L(1,B)P(B), if umbrella is taken 
= (loss if rain)(probability of rain) + (loss if no rain)(probability of) 
with unbrella with umbrella no rain 
= (0,50)(0.4) + (0.50)(0.6) 
= $0.50 
L(2) = L(2,A)P(A) + L(2,B)P(B), if umbrella is not taken 
= ( loss if rain )(probability of rain) + (loss if no rain)(probability of) 
without umbrella without umbrella no rain 
= (2.00)(0.40) + (0)(0.60) 
= $0.80 
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The action with minimum expected loss is action 1, take umbrella, and 
it is this action that the gentleman should accept. 


The logic for this strategy of choosing that particular action with 
minimum expected loss is that if this minimum expected loss selection 
criteria is employed in a large number of action selection problems, the 
actual total loss will be less than if any other method of action selec- 
tion were to be employed. 


Clearly, the gentleman has no mental anguish in making decisions 
once that he has decided to let the minimum expected loss strategy 
dominate his life. All that he requires is the probabilities of actual 
future conditions P(A), P(B), and the values of loss due to actions and 
Samoieaons®L(15A), L(2,A), L(1,B), L(2,B). 1/ 


A very interesting extension of the decision problem involves the 
selection of investigations which can reduce the uncertainty of future 
or unknown conditions. For example, is it worthwhile for the gentleman 
to invest $0.10 for a telephone call to the weather bureau, if their 
prediction of "rain" or "no rain" has a chance of 0.90 of being correct? 
Problems such as this are termed as "Pre-Posterior Analysis" since they 
involve selection of investigations previous to the "Posterior Analysis" 
of action selection, and the method of solution will be developed in 
Part II of this report. 


Now, however, civil engineering examples will be presented in order 
to extend the minimum expected loss concept to a more important class 
of practical problems. 


In all of the examples, the value or utility of actions and condi- 
tions will be expressed in terms of engineering estimates of cost or 
losses, since most engineers can easily understand how these values may 
be obtained. However, it should be understood that the minimum expected 
loss, or maximum expected gain concept can employ “utility values" which 
can represent social and environmental benefits on a rating scale which 
may differ from conventional monetary values, 


A Dam Spillway Capacity Selection Problem 


An earth dam is to be designed for a 100 year life. The designer 
must decide whether to select one of the following actions: 


1 
2 


construct 1000 cfs spillway capacity or 
construct 1500 cfs spillway capacity, at a cost of .5 million 
dollars more than the 1000 cfs spillway capacity dam. 


T/ If the nature of the problem requires estimating the losses or 
probabilities of future conditions then careful estimates are 
essential 
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actual future conditions are: 


A = Flood flow greater than 1000 cfs. (greater flow causes wash 
out of dam and complete destruction) 

B = flood flow less than 1000 cfs. 

C = flood flow greater than 1500 cfs. 

D = flood flow less than 1500 cfs. 


hydrological records provide the condition probabilities: 


P(A) = 0.2 
P(B) = 0.8 
P(C) = 0.1 
P(D) = 0.9 


and economic analysis of construction and flood damage costs provide 
the losses: 


L(1,A) = 10 million dollars 
L(1,B) = 0 

L(2,C) = 10.5 million dollars 
L(2,D) = .5 million dollars 


The expected losses are: 


L(1) = 0.2 (10) + 0.8 (0) = 2 million 
L(2) = 0.1 (10.5) +.09(0.5) = 1.50 million 


Here the minimum expected loss action is #2 = "construct 1500 cfs 
capacity." 
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(c) A General Formulation of the Decision or Action Selection Problem 


In this section, it is desired to formalize the decision making pro- 
cedure involving the selection of one particular action from a set of two 
or more alternative actions. The entire procedure will be developed with 
the use of the decision tree model. 


Again, it is best to present the general symbols in terms of a simple 
example as follows: A structural engineer is involved in the preliminary 
design of a reinforced concrete building. His decision problem involves 
the selection of one of the following actions: 

a} - Construct a concrete slab floor on two-way beams 
7 a 
{ ; 
a> = Construct a flat slab concrete floor without beams 
The decision of whether to take action a} or a, depends on the 


performance (strength, deflection, cost) of the strtictural floor system. 
The possible future conditions are: 


S1 - Perfect Performance of two-way beam slab 


Sy - Adequate Performance (slight vibration under load) 
of two-way beam slab 


S3 - Inadequate Performance (excess deflections and 
{S.J} = - cracking) of two-way beam slab 


Sq - Perfect Performance of Flat Slab 
Se - Adequate Performance of Flat Slab 
S¢ - Inadequate Performance of Flat Slab 
Now, the action selection must be made under uncertainty concerning 
the actual future conditions {S;}. The engineer must be able to obtain 


or assign probability values p(Ss) for these future conditions. Based 
on experience with previously constructed floor systems we have 


P(S}) = 0.6 

(So )e =) 0.0 Note P(S,) + P(S5) + P(S3) = 1 
for two-way beam slab 

P(S3) = 0,1 

P(S;) 2 

P(S,) = 0,5 

P(S;) = 0.3 NOtesP cman CS tar (S Jae 
for flat'slab  ~ 

P(S-) ony 
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Loss is incurred by the engineer according to the action a, he takes and 
the condition S, which eventually ome For example, ifthe action is 
ay» and S occuts, the loss L(a,, ) is the consequence. The engineer 
Ast be able to estimate the ies eatie of each action-condition come 
binations. 


Again experience with past construction projects provides the loss 
values for our example: 


L(aq.S}) = $8,000 
(here a negative loss indicates a profit) 
L(ajz.S2) = ~—- 42,000 
(positive values are losses) 
L(aj, S;) = L(a7.S3) = +10,000 
L(ao,Sq) = «10,000 
L(ao,Ss) = +2,000 
L(a2.S¢) = +8,000 
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L(a,,S,)=—8000 


AY) 
© P(S,)=0.3 


(¢} 
FA, % L(a,.S2)=+10,000 
0, %o ( j» 3) tT ® 
2 e L(apS4)=- 10,000 


Ss=adequate L(aa S5)=t2000 
P(Ss)=0.3 


ON L(a5S¢)= +8000 
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With this tree diagram we can follow through the entire possible condition 
and loss consequence sequence for any selected action or decision. 


A decision is thus defined as the taking of a particular action 
aj from the set of all possible actions, {a;}. The results or outcome 
of the action a; is the condition S; which qs random and has the proba- 


bility of occurrence P(S;). After taking a;, the S.; occurs and loss 
ae Sz) is incurred, The selection of the best action aj proceeds as 
folTows. 


With the loss values L(a S;) and probabilities Bion) Tots 


79 
possible to compute the expected value of each action a;, The engineer 
will then take the action which has the minimum expectetl loss. The 
expected loss of an action a. is found by the multiplication of the 
individual probabilities of Condition P(S;) by the corresponding 
loss L(a;, S;) for the action and the condition and then summing 
these produc Se 


L(a,) = E Lay, Sj) + P(S,) 
for all S, 
(where xy means the sum of all products for each S.) 
for all S; 
For the floor slab example, 
(ay) = = L(aqz, Ss) + P(S;) 
Be L(ayg057) * P(S)) + Lay, So) * P(S5) + (a), S35) © P(S2) 
=» (-8) « (0.6) + (2) + (0.3) + (10) + (0.1). = -3.2 
(where the (-) indicates an expected profit) 
L(ao) = £ L(ao, S;) ° P(S;) 
= L(ap, Sq) * P(Sq) + L(ap, Sg) * P(Se) + L(ao, Se) * P(S¢) 
pa |d)eem (0,5) + (2) (003) +.(8) (052) "= s=2.8 


Here the minimum expected loss or maximum expected profit action is aye 
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Evaluation of the Probabilities P(S;) of Future or Unknown 


Conditions or Situations 


In the previous examples the values for the probabilities P(S;) 
were found by: 


(1) Weather Department Experienced Judgement of the chances of 
rain, for the umbrella decision. 


(2) An empirical plot of a historical record flood flow data on 
probability paper for the dam spillway decision. 


(3) Past performance experience of the engineer with concrete 
floor systems allowed the assignment of the chances of levels 
of performance by judgement, in the slab decision problem. 


In most engineering applications the probability values will be 
assigned, as in this latter case, by the judgement of an expert with 
substantial experience concerning the future unknown condition. 


However, in this next example we will see that occasionally, 
probabilities can be computed in terms of accepted technical defini- 
tions, such as: "permissible error must not be exceeded with 90% 
reliability." 
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An Example from Surveying 


A civil engineer is engaged in a Surveying job concerning the 
running of a one mile length of levels. His decision problem is to 
select the use of either first order work or second order work. 
Specifically, the engineer's client requires that the one-mile 
length of levels be run with an error of closure less or equal 
to 0.015 feet. The level circuit must be repeated if a greater 
error occurs. A time limit does not allow more than one repeat 
to be made. 


a) Probability of Error: The final error of closure X is a random 
variable. This error is made up of the sums of the individual errors 
made at each rod reading in the circuit. A random variable such as 
X, which is a sum of a large number of random rod reading errors, 

has a "Normal" or "bell-shaped" probability distribution. 
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This curve is centered on the zero algebraic average value of 
error, representing the fact that error can be either positive or 
negative with equal likelihood. The spread of the curve is repre- 
sented by the standard error o, and probabilities of errors within 
given limits are provided by areas under the curve. For example: 


68.26% 95.44% 


IS. 87% 
E 
atts @).<p< 0, =2 One (@)u) t20n 
Prob [-o S<E<+o] = 68.26% Prob [-20 <E< 20] = 95.44% 


In fact, Normal Probability Tables (see Reference (2), page 555) can 
be used to find the area under any specified multiple of the standard 
error o. 


(l0O- _/)% 
or 
(I-.Z) x 100% 


lf 2 in decimal 


BL / 2A 
-K .O- +K-O- 
</2 O <L/2 


Prob [-K., 72° <E< +K 7204 = (1 = a)% 
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Thus with Normal Tables and a given value of o we can find the error 
limits +Kys20 such that the probability of remaining within these 
limits is a given value of (1 - a)% Or, for a given value of 

(1 - a)% value we can evaluate the standard error o. This is 
applicable to our particular problem: 


First and second order level work are specified procedures 
(see Reference (3) pages 9 and 12) such that for a 1 mile length 
of levels, 
First Order Permissible Error = Ej 
E} = 0.012 ft. 
Prob [-0.012<E<+0.012] = 90% 
This means that from experience, the error has been found to remain 
within +0.012 ft. within 90% of the time when First Order work procedures 
are employed. 
second Order Permissible Error = Eo 
E> = 0.025 ft. 
Prob [-0.025<E< +0.025] = 90% 
The above specifications for the (1 = a)% = 90% reliability of permis- 
sible error allows the determination of the standard errors o, and o, for 
first and second order work respectively. Normal Probability Tables give 
the multiple Ka/2 = K5y = 1.64 such that 

Prob [-1.640< Ex+1.640] = 90% 

and therefore for: 
First Order Work 


1.640, = E} = 0.012 


1 
giving the first order standard error, 


oA sae 
3, sy yf tla 0.0073 


Second Order Work 


1.6405 = Ep = 0.025 
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giving the second order standard error, 


Se Ue 
oe ACT a 0.0152 


With these o values the probabilities 


P[-0.01I5<E<+0.015] 


i 


can be evaluated from Normal Tables for first and second order work as 


follows, 


a/2 = .02 = 2% (From Tables) 
a -.04 
0.96 = 96% 


]ea 


OI 
Second Order 
0.015 0,015 | 
Ki/2 = 02 ~ 0.0152 = 9-99 
a/2 = 0.161 = 16.1% (From Tables) 
Oesehesee 
lea = 0.678 = 67.8% 


Defining the actual future consitions {S.}, (|E] is absolute value of E) 


le] S 0.015, 


lel S 0.015, 


le] >0.015, 


lel >0.015, 


with first order work, 


with second order work 


with first order work, 


with second order work 
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Defining the Actions {aj} 


First Order Work 


= 
a, = Second Order Work 
The Probabilities are as found above, 


0.678 


P(S,) =..0.96 P(S3) 
P(S.) = 0.04 P(S,) SEU oce 
b) The Loss Calculations: 


Both orders require 5 man crew: 


2 men at $5/hr = 10 

3 men at $5/hr = 9 

Wild Level at $l/hr = ] 
“$20/hr 


First order requires shorter shot lengths, i.e., more time 
time = 8 hrs at 20 = $160 


Second order 


time = 6 hrs at 20 $120 

Engineer appraises this job at $200 

Note: It is assumed that once the engineer begins with a given order 

then he will use that same order again if the first outcome is unsuccessful 
(E >.015) 

The actions, outcomes and losses are as follows: 


2nd outcome 


1f Ist 
Action Outcome unsuccessful Loss 
a S, (success ) L(a,5S,) = =200 + 160 = ~$40 (a profit) 
a, S, (failure) S, (success) L(a, S595.) = -200 +160 + 160 = $120 (a loss) 
a, S, (failure) S, (failure) L(a, 45,45.) = 160 + 160 = $320 
a, S.,(success ) -- L(a,4S,) = =-200 + 120 = -$80 
a, S, (failure) S,(success ) L(8954 9S) = «200 + 120 + 120 = $40 
ay S, (failure) S, (failure) L(4 6S, 9S4) = 120 + 120 = $240 
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& CO) 
& P(S o)= 04 


Sz = SUCCESS 
P(Sz)= .678 
S 


L(a),S, = -$40 


_ -(a),S5,Sa)° $120 


My 
oe ox Lia, 199299 )=$320 


Lta5, S3)-$80 


L(a,Sq,S3)*$40 
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d) The Expected Loss: 


L(ay) = (.96)(-40) + (.04) [(.96)(120) + (,04)(320)] = -$33.28 (an expected 
(success ) (failure) (success) (failure) profit 


L(a,) = (.678)(-80) + (.322)[(.678)(40) + (.322)(240)] = -$20.62 
(success ) (failure) (success) (failure) 


Choose a,, First Order Work, with minimum expected loss or maximum 
expected profit. 
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TI. DECISIONS CONCERNING THE SELECTION OF AN INVESTIGATION 
| (PRE-=POSTERIOR ANALYSIS) 


We return now to the umbrella problem given in part I, by asking whether 
it is worthwhile to invest $0.10 for a telephone call to the weather 
bureau, if the probability of their predicting rain when indeed it 

will rain and the probability of their predicting no rain when indeed it 
will not rain are both .90. 


Before solving this problem we will introduce some needed notation 
(1)-(4) and formulas (5)-(7). 


(1) P(AB) = The probability that both A and B will occur. 
(2) P(A) = The probability that A will not occur. 


(3) P(AB) = The probability that A will occur and B will not occur, 


(4) P(BJA) = The probability (conditional probability) that B will 
Occur given that A has occurred, 


P(B|A) may be found by determining P(BA) and P(A) and then dividing the 
former by the latter as shown in formula (5) 


(5) P(BJA) = r(BA) 


Given P'(A), P(B) and the conditional probability P(B|A) the "reverse" 
conditional probability P"(A|B) may be found by formula (6) below: 
(See reference (1) p. 119) 


(6) P"(A|B) = P(BJA) P'(A) , where superscripts (', ") are used to 
P(B 
illustrate the following differences in these probabilities. P'(A) 
(usually Stains is the original (prior) probability that A will occur, 
whereas P"(A|B) is the revised (posterior) conditional probability that 
A will occur given that more information has been gathered and it has 
been found that B has occurred. 


Finally formula (7) below states that the probability that A will occur 
is the sum of the probability that A and B will occur and the probability 
that A and. not B will occur. 


(7) P(A) = P(AB) +  P{AB) 


We have seen that if no caTl has been made to the weather bureau then the 
best course of action, based upon the minimum expected loss ($0.50), is 
"take umbrella". We now determine what the expected loss will be if a 

call is made to the weather bureau followed by the best further action after 
their report. We then choose to call the weather bureau if this expected 
loss is less then $0.50. 
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Below are tabulated the given probabilities. 


Prior Probabilities 


P'(R) = 4; the prior probability of rain 
P'(R) = .63 the prior probability of no rain 


Conditional probabilities we wish to "reverse" to obtain P"(rain|report rain) etc. 


P(RR|R) = .93 the probability of reporting it will rain when it will. 

P(RR IR) = .1; the probability of reporting it will not rain when it will. 
P(RR|R) = .1; the probability of reporting it will rain when it will not. 
P(RR|R) = .9; the probability of reporting it will not rain when it will not. 


Below we determine the unknown probabilities. 


The following two probabilities needed for determining the "reversed" posterior 
probabilities are: 


P(RR); the probability of reporting it will rain and 


P(RR) the probability of reporting it will not rain. 


They are as follows: 


P(RR) = P(RReR) + P(RReR) 


= P(RRIR)P'(R) + P(RRIR)P*(R) 
= (.9)(.4) +  (.1)(.6) 
= 42 

Hence (RR) = 1-P(RR) = .58 


We now determine the posterior probabilities of rain or no rain given 
the weather report as follows: 


P"(R|RR) = the probability of rain given a report of rain 


P(RRI|R)P'(R) _ (.9)(.4) . 
St ee 857 
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Hence: 


P" (R|RR) | 


P"(R|RR) 


Hence: 


P"(RIRR) 


the probability of no rain given a report of rain 
1seP"(RIRR) = .143 


the probability of rain given a report of no rain 


P(RR|R)P*(R) (.1)(.4) 


= = = .069 
P(RR) (.58) 


the probability of no rain given a report of no rain 


1 = P"(RIRR) = .931 
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L(a,)=L(a),1)= $0.50 


2 > ly 
AS L(a,,,RR,LR)= $0.60 
eK 
< t 
ON Ee est “ee 
ig gS 
we f 
F) ye Pa 2 
% 5 v AQ L(a,,RR,1,R)=$ 0.60 
rs % of oe) 5 a 
5 oo *y 
%\ Oe” Fey "ea, 
slo" OS a 
SSS 4 
"oe xe wendrers u RR,2,R)=$2.10 
Rise Vin Lan \= 4a L(a5, 9G» ° 
} AO 
*- Se Seg ae “oe 
Up Dog ee 
a, Lae 
be 20 re L(a5,RR,2,R)=$0.10 
? E>” 
-), 
‘S93, 


V/ L(A5. RR, 1, R) is the loss resulting from calling the weather bureau 


and obtaining a report of rain and taking the unbrella when it indeed 
rains. 


2/ Note: These two actions are those ruled out by using the principle 
of least expected loss at that action juncture. You may compute these 
two expected losses as an exercise. The significance of this is that 
for each possible result of every investigation candidate, the best 
choice of further action is both (1) necessary for selecting the best 
method of investigation and (2) available for immediate use after the 


results of the best investigation. This final step is called terminal 
analysis and is discussed later in this section. 
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The expected loss L(a,) resulting from investigation ay "Call the weather 
bureau" is: 


L(ay) = P(RR)EP*(R|RR)L(ap4RRyT4R) + P"(R[RR)L(ay4RR, TR) 


+ 


P(RR)[P"(R|RR)L(apyRRy2,R) + P"(R|RR)L(agyRR42,R) ] 


(.42)[(.857 x $0.60 + (.143) x $0.60)] 
+ 

(.58)[(.069) x $2.10 + (.931) x $0.10] 
$0.39 


Hence we conclude that it is best to call the weather bureau based 
upon the principle of least expected cost. 


AN EXAMPLE INVOLVING CHOICE OF A GEOLOGICAL SITE INVESTIGATION 


In the preparation of a road construction contract proposal, the 
engineer must decide which type of cut slope to Specify for a 1/4 

mile length of hill-side cut. The type of cut slope depends upon 

the unknown orientation condition of geological strata at the cut; 

where differing orientations of strata can cause different rock slide 
problems on the cut slope. The engineer can obtain information about 
the unknown strata by ordering one of two levels of geological site 
investigations, and therefore the pre=-posterior decision analysis is 

to select the optimum type of site investigation. Detailed descriptions 
of the actions, orientation conditions, investigations, results, losses, 
and probabilities are given as follows. 


A. Set of Actions {aj} = {a} ao} 
aj = 1:1 Cut Slope (to be called for in the contract proposal) 


A5Aa inci: 1 Cut Slope 


Cut volumes 

are substantially 
different for each 
type of specified 
Slope. 
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B. Orientation Conditions tog £51;.-Soi 


S} = Orthogonal (to cut) Strata 


Note: Can use aj = 1:1 
slope for good 
slide control 


So = Horizontal Strata 


Note: Can use a> = 2] 
slope for good 
Slide control 


C. List of Possible Site Investigations {e,} = {@,5 eo} 


Bond Hire geological consultant with manual sub-surface exploration 
by shovels and augers. 


Boy 2 Hire geological consultant with drilling rig for deep cores. 
D. Set of Possible Results of Investigations {Z,} SietZ7 scot 

Z, = Indication that S, exists 

Za = Indication that S» exists 


Note: An indication does not definitely establish that a given condition 
exists. It may, however, increase or decrease the probability of 
the presence of a given orientation condition S,. 
E, Probability Assignments 
1. Prior Probabilities for the Unknown Orientation Conditions: 


In order to select one of the two types of investigations e,, the 
engineer must have some idea of the probability of the existence of 


26 


org 
4 


‘ 


(2 al = (2h anmtetbeed aolsasi 


‘ri tdae 


ytyob 


ey] 


‘~ 


al | 


OS f 


sian? (dun edd Tenogotd 


SE TEL id 2aveat sjfe2 aternnn Yo 


w Feat y2nes fea botele avih * 
esoove bas 2eTevat2 ¥d 


Ww jcetipenog Testeefose ari = 
utteovnl to eglugoh atdteen4 te 
2tetxs 2 todd hobgeolial #9 
2iehas ef jetd netygentinl = 


Pal 


yis?tahtsb J60 each watz aatoat 
"9 S2soront ,TevewoR , vam FF oe 
'setnsive nevib #40 Soneaeag ang 
zane teh ei 


rwondcl) say 10% abs URaaa vor. 
eo? ows. wi To one ae 
eq seb Yo aah semble toa ens 


os 


each unknown strata condition Sj. Now, since the engineer assigns a 
number value P'(S;) for these probabilities prior to the investiga- 
tion, they are termed as "Prior Probabilities of Conditions" = P'(S;). 
Here it is necessary for the engineer to employ all of his technical 
knowledge and experience to best evaluate the prior chances P'(S;) 
before having any knowledge of investigation results. Of course, 
there will be cases where the engineer or decision maker will be 
totally ignorant of the conditions. This state of total ignorance 
provides a base line for probability assignment; for example if the 
engineer has no knowledge about orientation conditions S] or So, he 
should assign equal probabilities, totaling up to unity, 


BCs Mea cere So) = = al {2 


That is, total ignorance may be expressed by equal chances for all 
unknown conditions. 


In our problem however, the engineer has seen the character of the 
terrain in his past travels and has a slight preference for the 
horizontal orientation S9, such that his prior probability assign- 
ments are P'(S,) = 0.4, and P'(So) = 0.6. 


2. Conditional Probabilities of Investigation Results for a Given 
Investigation and Actual Condition: 


In many types of engineering investigations, technical knowledge and 
experience will allow the assignment of a Conditional Probability 
P(zy|S;, e,) of obtaining a certain result z, from a given investi- 
gation e,when an actual known condition So exists. The term "condi- 
tional probability" implies that the probability is valid only for the 
condition that S; is the true existing Situation. 


The given condition S; allows us either to draw physically, or des- 
cribe mathematically, the true given condition S;, and then estimate 
or calculate the chance of obtaining a given result z, with investi- 
gation e,. For example, the geologist is able to draw the actual S 
condition of horizontal strata on the roadway cut cross-sections, 
then show the types of his investigations e, or eo; and then estimate 
the chances of getting z} or Z5 from each aE these investigations; 
where natural variability of the horizontal strata might yield er- 
roneous results due to just "bad luck" at the investigation sites. 


Hand Dug Ditch Investigation e, may yield 


erroneous 2 indication here due to coinci- 
dence of upward strata variation at the 


ditch location. Chance of this "bad luck" 
coincidence is estimated as 


P(z)|Sos ey) = 0.3 
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Drill Hole Boring Investigation eo may 
yield erroneous 24 indication here. Chance 
of this "bad luck" coincidence of a boring 
with an upward variability of strata is 
estimated by the geologist as 


P(z,|So. Oo) Pcie Gel 


The geologist is able to assign this estimated probability value because 
he has the knowledge of the frequency of “upward strata irregularities" 
that may occur when the true strata is horizontal So. 


For the assignment of P(z [Sis er)» in general, the investigator asks 
himself the following question: “What is the probability that my in- 
vestigation ey will yield a result z, when the true condition is S.;?" 
or "With what reliability will my in¥estigation e, detect the condition 
S; when it truly exists?" The fact that these conditional probabilities 
can be asSigned by an experienced expert, with access to past records 
and technical knowledge relating to the given actual conditions and 
investigations, is a most important input phase of the computations 
leading to the pre-posterior selection of an optimum investigation. 
This is because the conditional probabilities R(zs) Sis ey) represent 

a numerical measure of the investigator's opinion Concerning the ac- 
curacy or quality of the results of a given proposed investigation e,. 
For example, if the experience and knowledge concerning a certain in- 
vestigation e, is such that the investigator feels very confident that 
it will furnish a reliable indication z, of the condition S. when S. 
really exists, then he may express this confidence numerically by as- 
signing a high value for P(z [Sis e,) such as 0.9. In the analysis 
procedure, this high P(z-|Siy e.) Value will favor the actual selec- 
tion of e, as the optimum investigation, if its cost is not too high. 
The complete set of Conditional Probabilities are given below: 


method of investigation 


indicated 
condition 


t 


actual “condition 


For example: the geologist knows by experience that there is a 


0.7 chance of zj when investigation e, is made and Sj is the 
actual condition. 
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Losses or Costs L (ek, Zo, aj, Sj) for each Combination of Investiga- 
tion, Result, Action, and Neeual Condition. 


All costs are expressed in terms of Equivalent Annual Cost in Dollars. 


Construction Costs: These vary for each action a; because of the 
ifferent earthwork cut volumes for each type of dpecified Slope. 
These costs can therefore be estimated from known volumes and unit 

COS tS. 


Construction Cost for aj = CC (aj) 


CC (a7) = 0, here the cost for this minimum necessary action is 
used as the zero base. 


CC (a9) = 200 


Maintenance Costs: These vary for each combination of specified slope 


a+, and actual orientation condition Sj. Here, particularly, past 
available maintenance records for each of the above combinations, can 
provide cost data. 

Maintenance Cost for aj, with Sj = MC (ay, S;) 

MC (a), Sy) = Q; MC (ao, S}) = 0 


Because each slope can adequately control slides with the S] orienta- 
tion condition. 


MC (ay> So) = 400 ; MC (ao, So) = 0; 


Because hee 1:1 slope does not control slides with the So Orientation 
condition. 


Investigation Costs: These vary with the manpower, equipment and time 
required for each type of investigation e,. Here it is assumed that 
the cost of any investigation e, does not vary with the result or 
indication Zz, detected by the investigation. 


Investigation Cost for investigation e, with result z, = IC (e,. Ze) 


IG (ey, Z) = 100; IC (Qo, Zo) = 200 
All are valid for either result Z1» Zoe 


Total Cost or Losses: The total loss for a given investigation eps 
result z,, action ays and condition S; is 


ahs 


L(e,s Zoo ass S;) = IC(e,, zo) + MC(ai. S;) + CC (a, 
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This total cost or loss is given in the following table, 


Ller,s Zos ays S;) ay ao 
L(e}» Zo» aj» So) 500 300 
L(eo, Zoe aj» S}) 200 400 


The above costs do not vary with investigation results 215 Zo. 
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A Flow Chart Guide for the Reader 


The following Flow Chart is provided to Serve as a general 
guide to the reader, during the simple, but lengthy, calculations 
involved in this numerical example. Definitions of all symbols 
will be given directly in the example. 


It should be recalled that the basic objective is to decide 
which type of investigation e,, should be made in order to improve 
the unknown condition probability values from the rough prior (prior 
to the investigation e;,) values of P'(S;) to refined or sharper 
posterior (posterior to having conducted e, with a result of Zp) 
values of P"(S;|z,, e,). For example, the engineer wants to know 
ieeit is worthwhi Ye to spend $200 on the boring investigations e9, 
which has the reliability P(z]|S}, e2) = 0.9 of providing an 
indication z} of upward strata when upward strata S} is the true 
condition; such that the prior probability of upward strata 
P'(S}) = 0.4 would be changed to the sharper posterior value 
(to be calculated in the example) of P"(S]|z]» e2) = 0.86. 
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FLOW CHART 


FOR PRE=POSTERIOR ANALYSIS 


Given the Input Quantities 
Basa} 


P(Z. |Sq sex) 


L( ek 929 92; 94) 


Compute Probabilities 


P(z, le.) = Pe OP( Zoi Sy eer) eeehaiS4) 
Sum 
for all S; 


P(z,[Ss,e,) . P(S 
PM(S;|Z,2e,) = Ze [ek 


Compute Expected Losses for Actions 
L(e, 52, 4a,) = uy L(G, sZ prays S;) : P"(S.[Z,se,) 
for all S; 


then select the minimum value 


* = ° ° ae 
L(e, 4Z, 9a ) = minimum of L(e, Za.) 


Compute Expected Loss or Cost of Investigation 


L(e,) = E L(epsZpsa*) P(Zp ley) 
sum 


for all Zo 


Then select the minimum value 


L(e*) = minimum of L(e, ) 


The optimum or minimum expected loss investigation 
is that e, which has L(e*) 
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Calculation of Further Required Probabilities 


In order to proceed with the selection of the optimum, "minimum 
average loss", investigation, it is necessary to calculate: 


(1) The Probability P(z,|e,.) of obtaining a result zy from a given 
investigation e,. 


(2) The Posterior Probability P'(Si|Zp. e.) of the existence of 
the true condition S; when a given investigation e, provides 
a given result Zp. 


Both of the above probabilities can be calculated from the available 
"judgment assigned" probability values. 


P'(S;) = Prior probability of a true condition Sj. 

P(zp|Sis ex) = Probability of obtaining a result z», from a given 
investigation e, when a given true condition S; 
exists, 


The required probability calculation procedures are best presented 
by a Venn diagram. (See Reference (1), page 14) 
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The joint probability P(S;, zpg|e,) of the occurrence of both the true 
condition Sj and the result Zz, of a given investigation e, may be 
found by use of the basic definition of conditional probability, 

(See Reference (1), page 62), 


P(S. » Z, |e) 


BRR aL Uaaeialianys Try 
which gives 
= A ] 
P(S., z, le) = P(z, |S. e) p (S,) 


in terms of the available probability values. 


From the Venn diagram the inner area which represents P(zgle,) is seen 
to be equal to the sum of all the joint probabilities, 


2 


thus, uSing the above equation for P(S., z les 
j 


| Wats Jet aS) 
J 


P(z le.) = P(z |S. e.) sahue losis 


j 1 


Now if we again employ the definition of conditional probability, 
the Posterior probability is given by 


P(Ss» Z,/e,) 
P'(S.|z,, e,) 3 P(z le) 


using the above results, 


P(z,|S;> e,) 3 P*(Ss) 


P (S.|z05 e.) 


fl 1 DO 


P(z |S. e.) . P'(S.) 


which is in the form of Bayes' Theorem. (See Reference (1), page 119). 
For our example, above probability computations are organized as follows: 
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Selection of the Investigation with Minimum Average Loss by 


Working Backwards on the Decision Tree 


The decision tree display of the investigation selection problem is: 
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S\ 


L(@),Z25a9,S)) 


L(e |? 259908 : ) 


een oneion e, has its own respective set of result probabilities 
P(z,je 
mi Ks e 


The result z,|e, improves the knowledge concerning the true but unknown 
condition Sj, this improvement in knowledge being from the prior proba- 
bility P'(s;) to the posterior probability P"(S4|zp, ek), which is used 
to select the decisions aj. 


Now since, before the selection of any investigation ex, all information 
is known only on a probabilistic basis such as P(zglek) and P"(S4|Z1, ek)s 
the value of an investigation e;, must be computed on an expected or aver- 
age cost basis. Decisions concerning the selection of actions and inves- 
tigations will be based on the minimum average loss or cost basis. This 
computation proceeds as follows, by working backwards from the action 
selection branch of the decision tree: 


L(e,..Z . 10; 8; 
P sop ize k) 
Position of Engineer Position of Engineer when 
when computing L(e,) computing L (ex, Zg, aj) and 
and finding minimum finding the a* which provides 
L(e*) the minimum value L(ek, Zp, a*). 
Here, he imagines that z,|e, 


is known. 


For each of the combinations of investigations and corresponding results 
zglex, the engineer imagines that he has this information, and with the 
resulting posterior probability P"(S;|Zps ex)» he positions himself at 
the action selection branch point of the decision tree. At this point, 
he is able to compute the average loss of each action aje 
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For each of the combinations of investigations and corresponding results 
Zeleks the engineer imagines that he has this information, and with the 
resulting posterior probability P"(S;|z), e.), he positions himself at 

the action selection branch point of the decision tree. At this point, 


he is able to compute the average loss of each action aje 


aot 


L(evs Zo. a) = ; 


Wt 3 


A iT} 
Lleys Zp As S.) p (S. |Z). e,) 


and select the action a* with the minimum average loss L(e, Zea) + 
Thus for each of the possible investigation results z ley» thefe is a 
corresponding minimum average loss value L{e,, Zoe a*7, 


The engineer now moves to the investigation selection branch point of the 
decision tree and imagines that he has selected a given investigation e 
but has not yet obtained the results z_|e,. He does, however, know thekr 
probabilities P(z le). and the average Ibss or cost of e, can be computed 
by 


L(e,) = z L(er, Zuw a") P(z, ley) 
Q= 


This may be interpreted as being the average cost of all the optimum 
actions a* having minimum average loss, when investigation e, is used. 
The engineer, should select the investigation e* which has the minimum 
cost value L(e*). For our example, the above average loss computations 
are organized as follows: | 


om 


A. Compute L(e , z , a.) and find minimum = L(e ,z , a*) 
ocr 3 ky ad 
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Action aj 


u 
Pp (S., eps z,) 


2 
22 86 07 


a 
P"(S,) Bi 


e e e 


L(S,) 
L(S,) . P"(S,) 
P"(S,) 
L(S,) 


L(S,) . P"(S,) 


oy = 


aN peace 
F 


L(e, 522587) 


L(e, sZ, .a*) = 


Pritt ae 256 --- 
min L(e, 5Z)a5) 256 
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Action ao 


PN(Silexs 29) ey e2 
L(e@, 92, 28595; ) zy 25 zy Zo 
p*(s,) .61 a2? . 86 .07 
L(S,) 300 300 400 400 
L(S,) + p(s.) 183 66 344 28 
p™(S..) .39 .78 14 93 
L(S,) 300 300 400 400 
L(s,) > P*(S,) 177 234 56 372 
ee 
over all oF 
Be {e,s Z,s ap) 300* 300 400 400 


Le, zs a*) = 


ink = 300 — 400 
min L(e, Zs a) 


*These values are uniform because action ao provides adequate slide 
control for both conditions S, and So. 
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B. Compute L(e,) and find L(e*) = min. of L(e)s 


= 2 = Z 
where L =f *) oP Bais CURE EAL 
(e,) Pauly scrsa ) (z le) a (a*) P(z ) 

P(z 

1!*) e e 
L(e,.5Z,.a*) g 
P(z, |e,) 460.42 
L(e, 2, a*) 256 256 
L(a*) «P(z,) 118 108 
P(z, le, ) «04 258 
Le, 2,3”) 300 400 
L(a*)+P(z.,) 162 230 
2 —<_ 

z L(a*)+P(z,) 280 340 
g=] 


L(e*) = min L(e,) = 280 = L(e,) 


Therefore select e., the site investigation by the manual ditching 
method. | 
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C. 


POSTERIOR OR TERMINAL ANALYSIS, AFTER THE INVESTIGATION 


Once having selected the optimum investigation e*, the engineer orders 
its performance, and obtains the result zp. With this new information, he 
revises his prior probabilities P'(S;) to the posterior values PAS Het ak 
and proceeds to select the best action a; based on minimum expected loss 
calculations which employ the new revised P"(S|z,,e*) values. It is use- 
ful to note that all necessary calculations for this selection of the best 
action have already been performed during the pre-posterior analysis select- 
ion of e*; the expected loss Clez,s Zes aj) based on the corresponding revised 
"(Si|Zp, e,) values for all possible combinations of Ck» Ze, aj are all 
avai lab e in the calculation tables. 


For example, in terms of the roadway slope selection problem: 


The engineer, having selected e], now performs his investigation and 
obtains the outcome zp, here for example: Zp = Z}, an indication of the 
upward strata_condition S}. He will select the action with the minimum 
average loss L(e}, Z}, a*| which is found in the computation tables as 
Action aj with an average loss value of 256, Therefore select action aj, 
the 1: 1 slope. 


The entire decision sequence, as illustrated in this example, of select- 
ing the investigation e;, then obtaining the result zj, and finally selecting 
the action aj, is the best strategy that the engineer can employ in the face 
of uncertain conditions, It should be realized that his selection of action 
a}jhas a chance of being wrong, since there still remains a probability 
P"(Solez, 21) = 0.39 of having the horizontal strata condition So in the 
selected 1 : 1 slope of action aj, with the corresponding slide problems and 
high loss L(eq, Zs a}. S2) = $500. 


However, without knowing the true strata condition S;, the engineer has 
minimized his expected loss when he chooses a}, and if he employs this 
strategy of minimum expected loss in all of his many professional decisions, 
his total actual loss will be smaller than if he employed some other strategy 
or no decision making strategy at all. Further developments and examples 
may be found in references (4) and (5). 
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III. 


AN ADDITIONAL EXAMPLE 
Terminal Analysis (Decision analysis after the sampling investigation) 


Terminal analysis refers to the posterior analysis of the decision tree 
when new information about the conditions has been obtained by a sampling inves- 
tigation. The optimum action is selected after the sampling investigation has 
been performed. In this type of analysis, the cost of sampling is considered to 
be an already incurred cost (sunk cost) and no longer enters into the problem. 


The sampling investigation provides additional information concerning the 
unknown conditions and terminal analysis treats the problem of revising prior 
probabilities as a result of this sample information. The expected values of the 
various actions are computed using these revised probabilities (termed as poster- 
ior probabilities) to determine the best action. 


The example in this section will deal with a normal probability distribution 
which is found to fit a number of common conditions in engineering decisions. 
Statement: The reader will perhaps not fully understand some of the terminology 

and procedures -- however, they are standard methods of engineering statistics. 
This example is given to show the powerful capabilities of the method to solve 
practical problems ; lack of some procedural comprehension can easily be corrected 
by further study in standard texts. The normal distribution is fully tabulated 
in terms of the mean X (average value of the variable X) and a measure of dis- 
persion or scatter of the individual values of the variable around the average 
value known as the standard deviation o. The normal distribution has a prob-« 
ability density function which is a smooth, symmetric, continuous, bell-shaped 
curve as pictured in the figure below. The area under the curve over any inter- 
val on the horizontal axis represents the probability of the random variable, X, 
taking on a value in that interval. The area under the total curve is equal to 
unity. This curve reaches a maximum at the mean of the distribution. One half 
of the area lies on either side of the mean. The standard deviation o controls 
the spread of the curve; for example, a large o means a very wide curve. With 
any normal distribution, 0.50 of the area lies within + 0.67 standard deviations 
from the mean; 0.68 of the area lies within + 1.0 standard deviations; and 0.95 
of the area lies within + 1.96 standard deviations. 
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Example: An engineer is faced with the problem of accepting (action a]) or reject- 
ing Rectan a>) a particular batch of Class "A" mix concrete delivered to an im- 
portant job. The desired or specified mean concrete strength is 4000 psi. Assume 
that extensive observations on strengths of concrete of this Class "A" type have 
Shown that the "population" strength X is normally distributed with a standard 
deviation of oy = 900 psi. The "population" here refers to the whole class of 
concrete about which conclusions are to be drawn. The engineer, based on his 
experience in concrete mix design, assigns the following prior probabilities of 
obtaining significantly different actual mean values Sj when class "A" mix is 


specified. 
Actual Mean, Sj psi Prior Probability P'(S;) 
3500 0.1 
4000 (The desired 0.6 
value) 
4500 0.3 


Later, data from a sampling investigation will be used to revise or sharpen 
these prior probabilities. The data will consist of 9 cylinder tests with 
an arithmetic average value of 


An 3700 psi 
this observed average is an estimate of the true mean S.. 


First however, we shall perform a "prior" decision analysis with only the 
prior probabilities P'(S;). 


The engineer has studied the economic value problem extensively and has 
arrived at the following dollar losses for each of the conditions represented 
by the various possible actual means: 
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Actual Mean S; Receee Reject 
3500 $1000 $ 60 
4000 $ 0 $ 500 
4500 $ 0 $1000 


Using the prior probabilities, the expected loss can be computed as 
shown below. 


Prior 
Loss Expected Loss 
Actual Mean Sj Accept Reject P*(S;) Accept Reject 
3500 $1000 or $ O x Q.1 = $100 $ 0 
4000 $ 0 or $500 x 0.6 = $ 0 $300 
4500 $ OFT or ST000 xe OCSe = Se e0 $300 
Total expected loss $100 $600 


It can be seen that if the sample data is not used, the prior expected losses of 
accepting, or rejecting are $100 and $600, respectively. Thus, using prior prob- 
abilities only, the mix should be accepted. It is assumed that concrete cylinder 
strengths X follow a normal distribution; the parameters which "sufficiently" des-= 
cribe this process are the sample mean X and the number of samples n. For large ~ 
samples the distribution of sample means X is also normal and the standard devia- 
tion of the sample means from the true mean is given by 


where 


standard deviation of the cylinder strength X 


number of observations of X = size of sample 


— 
u 
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Normal Distribution of Cylinder Strength 
about a possible actual mean value S, 


Normal Distribution of a Sample Mean X 
equal to the arithmetic average of n cylinder tests X 
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In this example: 


poe = 300 psi 


X A/9 
The sample likelihood* or the probability of obtaining the observed sample mean X 
as a function of the true mean S; is computed for each possible value of the true 


X- S, 
Mis calculated, then using this value, the ordinate of 


means Sie First 
Co@e 


the normal curve is read from the tables of the normal distribution. (a sample 


is attached). 


inflection point 


standard 
Te 


true 
Ox ordinate 
ordinate 


“AE 


Standard Variable Sample Mean X 


* See Benjamin and Cornell, "Statistics and Probability for Civil Engineers," 
McGraw-Hill, for terminology and notation used in the following calculations. 


49 


nee alume bewieeedo ut? onintezde to ysl idedotg si te *poodt feat i" 


vit ocd 76 culsy sfdteaon rides 10% betuqmm> ef pe Neem sung a“ 
- 
Yo ofanibuo off ,owlev z?dd prifeu wot? ,baisiiwoles 2 Ian 
at 
; x 
cs) ,naltudisdeth fearon sit Ya eafdsd Sao mot bese oP s 


a 


=e 


",zvawmipad Itvtd vot @ifidsdart bas eotseteed2* af = tmstagd. 
.2nohjelyalseo potwotiol af3 wt bszu nofgsson ones . a un 


Ordinates of the Normal Curve 


1 2° 
Ub eee oy 
on 
z 00 07 .09 


The above ordinates give the probability density for z = (x-u)/o deviations 
from the mean yp (that is, z = 0). To fit a normal frequency curve to observed 


data consisting of n observations, multiply the ordinate from the table for 
any value of z by n/o. To fit a normal probability curve multiply the ordinate 
by l/o. 


Table A-3 "Basic Statistical Methods for Engineers and Scientists," 
Neville and Kennedy, International Textbook Company. 
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Condition 


3500 


4000 


3700 = 4500 


4500 


Condition 


4000 


4500 


standard 
variable 


X- Sj 


(0) 


X 


3700 = 3500 


300 


3700 = 4000 


300 


300 


Standard 


ordinate 


y 
(from Tables) 


ieay) 0.3187 
-1.00 Oe 242 
=-2.6/7 0.0113 
Sample 
likelihood Joint Measure 
2(Xsn|S.) 2(Xons5, |e) 
0.3187 0.03187 
300 = 300 
0.242 y 0.14520 
300 = 300 
OAT = 0.00330 
300 "’ 300 
r= 0.18043 
300 
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Posterior Prob. 


p(S.|Xsnse) 
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The sample likelihoods computed above are the relative probabilities of 


finding the actual sample mean given that the mean is one of three possible 
values 3500, 4000, and 4500 psi. 


The posterior probabilities are calculated by the following expression: 


P"(S;|z) = P'(S3)e2(z]S;)en(z) , where z = sample information (Xsn) 


or 


Posterior Probability = Prior Prob. times sample likelihood times 
a normalizing constant 


The joint measures are calculated by multiplying the prior probabilities 

of state by the sample likelihood given that particular state. The joint mea- 
sures are then sunmed to give the reciprocal of the normalizing constant. The 
normalizing constant is introduced so that the sum of the posterior probabilities, 
for all of the conditions considered, will equal unity. The posterior proba- 
bilities then are calculated by taking the ratio of each of the joint measures 

to their sum. Terminal analysis concludes with the computation of the expected 
posterior losses as shown in the following table: 


Posterior Posterior 
Condition LOsSau Probability Expected Loss 
S; Accept Reject P"(S.|Xyn,e) Accept Reject 
3500 $1000 $0 0.177 $177 $0 
4000 $0 $500 0,805 $0 $402 
4500 $0 $1000 0.018 $0 $ 18 
Total expected posterior loss = $177 $420 
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$1000 


SN 4° 
PRIOR EXPECTED As 
Loss 100 ay S5* 4000 
at (S2)= 0.6 $0 
: As se 
2 ($ %S0, 
wv ROA (@) 
6 3 $0 
MINIMUM PRI Ov 
EXPECTED LOSS $100 $0 
Sa A Qe oe 
¥ acy Fo 
i <c> SN K> 
& ny CON 
‘VO PRIOR EXPECTED er Sinaia 
ins Loss $600 \_ P%S,) = 06 $500 
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MINIMUM POSTERIOR # ee ag00 pitt 
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yy S\ (S\ 
Ph: 
POSTERIOR Earecie od tala sae $500 
Loss $420 P"(S,)= 0.805 
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GROUP OF TEST RESULTS GROUP OF PROBABILITIES & LOSSES 
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It should be noted that the posterior probabilities differ from the prior 
probabilities by virtue of the sample results. The sample mean strength was 
lower than the desired mean strength; this produces an increase in the expected 
cost of acceptance of the mix. It should also be noted that the limitation of 
S; to three specific values restricts the posterior probabilities to those same 


values. 
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